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ABSTRACT 

Microlensing observations towards M31 are a powerful tool for the study of the dark matter popu- 
lation in the form of MACHOs both in the Galaxy and the M31 halos, a still unresolved issue, as well 
as for the analysis of the characteristics of the M31 luminous populations. In this work we present the 
second year results of our pixel lensing campaign carried out towards M31 using the 152 cm Cassini 
telescope in Loiano. We have established an automatic pipeline for the detection and the character- 
isation of microlensing variations. We have carried out a complete simulation of the experiment and 
evaluated the expected signal, including an analysis of the efficiency of our pipeline. As a result, we 
select 1-2 candidate microlensing events (according to different selection criteria). This output is in 
agreement with the expected rate of M31 self-lensing events. However, the statistics are still too low 
to draw definitive conclusions on MACHO lensing. 

Subject headings: dark matter — gravitational lensing — galaxies: halos — galaxies: individual (M31, 
NGC 224) — Galaxy: halo 



1. INTRODUCTION 

The search for microlensing events aimed at the char- 
acterisation of the MACH O distr i bution in galactic ha- 
los, first discussed bv iPaczvhskil (|1986l ). is by now an 
established technique. The results obtained up to now 
are, however, debated. Towards the LMC the MA- 
CHO group have claimed the detection of a MACHO 
signal from objects of ~ 0.4 M Q that would constitute a 
halo mass frac tion of about / ~ 0.2 ([Alcock et al.l 120001 
Bennett 2005), whereas the EROS group have found no 
candidate microlensing events and put a rather strin- 
gent upper limit on the same quantity, / < 0.1 in the 
MACHO mass range preferred by the MACHO results 
(jTisserand et alJ 12007( 1 . The issue of the nature of the 
detected candidates still remains an open question ( Sah 
1994 1 : |WullL99l iMancini et al.l[2 00l ICalchi Novati et al 
200a Evans fc Belokurovil200a 
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The contradictory results obtained towards the Magel- 
lanic Clouds challenge one to probe the MACHO dis- 
tribution along different lines of sight. Beyond the 
Galaxy, M31 represents the next most suitable tar- 
get for micro l ensing: searches (|Crottslll992l : iBaillon et alJ 
Il993t Uetzerl |l994). Looking at it from outside, we 
can globally study the M31 halo; the line of sight to- 
wards M31 allows one to probe the Galactic halo along 
a different direction; the inclination of the M31 disk 
is expected to give a clear signature in the spatial 
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distribution for microlensing events due to lenses in 
the M31 halo. Several obse rvational campaign s have 
been carried out: AGAPE ([Ansari et al.l [1997), who 
presented the first con vincing microlensin g candidate 
along thi s line of sight (| Ansari et al.l IT999T ) . Columbia- 
VATT (ICrotts fc Tomanevi I1996H, POINT- AGAPE 
(|Auriere et al.l 120011: iPaulin-Henriksson et al.l [2 003) , 
SLOTT- AGAPE (ICalchi Novati et all 120021120031). We- 
CAP P (iRlffeser et alJ 120031). M EGA (Ide Jong et all 
I2004D . NainiTal (|Joshi et al J 120051) . The detection of a 
few microlensing candidates has been reported, as well as 
first, though contradictory, conclusions on the MACHO 
content along this line of sight. The POINT-AGAPE 
group have reporte d evidence of a MACHO signal 
(|Calchi Novati et alJ [2005h . whereas the MEGA group 
have concluded that their detected signal is compatible 
with the expected M31 self-lensing r a te (Ide Jong et al.l 
I2006D . Very recently, iRiffeser et all (|2008|) have pre- 
sented a new analysis of a previoulsy reported bright 
event observed towards the M31 central region. Taking 
into account the effects of the source's finite size, they 
have concluded that the lens of this event should be at- 
tributed to the MACHO population. Finally, we recall 
that a few in teresting attempts have also been propose d 
(|Totanill2003T ). or already carried out (|Baltz et al.ll2004D . 
towards targets located beyond the Local Group. 

In 2006 we began a new observational microlensing 
campaign towards M31 using the Cassini 152 cm tele- 
scope at the "Osservatorio Astronomico di Bologna" 
(OAB) located in Loiano 8 . The res ults of the first year 
pilot s eason have been discussed in ICalchi Novati et ail 
((20071) . In this paper we discuss the second year cam- 
paign. As the main result, we have carried out a complete 
analysis of the microlensing flux variations, selected two 
microlensing candidates and compared them with the ex- 
pected microlensing signal. In Sect. [2] we present the 
observational setup and outline our data reduction and 

8 http:/ /www. bo. astro.it/loiano/index. htm. 
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analysis technique. In Sect. [3] we present our pipeline for 
the search for microlensing-like flux variations. In Sect.H] 
we present the simulation of the experiment with an eval- 
uation of the expected signal. In Sect. [5] we discuss the 
main results of the present analysis. Finally, in the Ap- 
pendix we describe in some detail some of the steps of 
our selection pipeline and of our Monte Carlo scheme. 

2. DATA ANALYSIS 

2.1. Observational setup, data acquisition and reduction 

The data have been collected at the 152 cm Cassini 
Telescope located in Loiano (Bologna, Italy). We 
make use of a CCD EEV of 1340 x 1300 pixels of 
0.58" for a total field of view of 13' x 12.6', 
with gain of 1.0e _ /ADU (this value has changed 
with respect to that of the first season because of 
some electronics problems) and low read-out noise 
(3.5 e _ /px). We have been monitoring two fields of 
view around the inner M31 region, centered respec- 
tively in RA=0h42m50s, DEC=41°23'57" ("North") and 
RA=0h42m50s, DEC=41°08'23" ("South") (J2000), so 
to leave out the innermost (~ 3') M31 bulge region, and 
with the CCD axes parallel to the south-north and east- 
west directions so to get the maximum field overlap with 
previous campaigns. This second year campaign lasted 
50 consecutive full nights, from November 11 to Decem- 
ber 31, 2007, with a fraction of good weather of almost 
60%. In order to test for achromaticity, data have been 
acquired in two bandpasses (similar to Cousins R and /) , 
with exposure times up to 6 minutes per frame. Over- 
all we collected about 410 (280) exposures per field over 
31 nights in the R (I) band 9 . Typical seeing values are 
~ 2" (somewhat worse than during the first season) . Sky 
flat frames were taken whenever possible so as to build a 
master flat image (per filter), and standard data reduc- 
tion, including bias subtraction, was carried out using the 
IRAF package 10 . We corrected / filter data for fringe ef- 
fects. The analysis presented in this paper is based on 
the 2007 season data only. 

2.2. Image analysis 

As for the preliminary image analysis we closely follow 
the strategy (th e "pixel-photometry") adopted by the 
AGA PE group (jAnsari et all Il997t ICalchi Novati et all 
2002), wherein each image is geometrically and photo- 
metrically aligned relative to a reference image. To ac- 
count for seeing variations we then substitute the flux of 
each pixel with that of the corresponding 5-pixel square 
"superpixel" centered on it (whose size is chosen so to 
cover most of the average seeing disk) and then apply 
an empirical, linear, correction in the flux, again cali- 
brating each image with respect to the reference image. 
The final expression for the flux error accounts both for 
the statistical error in the flux count and for the residual 
error linked to the seeing correction procedure. Finally, 
in order to increase the signal-to-noise ratio, we combine 
the images so to get 1 data point per night per filter. 

9 We do not have exactly the same number of data points per 
night per filter for the two fields, so that the indicated value is 
actually an upper limit. Furthermore, within the analysis, we ex- 
clude a small fraction of data points that show anomalously large 
relative error values, usually associated with poor seeing conditions 
or, more generally, poor image quality. 

10 http://iraf.noao.edu/. 



We evaluate the calibration zero point for the instru- 
mental magnitude versus standard (RI)c magnitudes by 
using a sample of secondary reference stars (|Massev et al.l 
I2006fi . We find for R and I bands data C R = 23.1 and 
Cj = 22.7, respectively (the reported values correspond- 
ing to the standard magnitude for an object with instru- 
mental magnitude of 1 ADU s _1 ). 

3. MICROLENSING EVENT SEARCH PIPELINE 

We have established a fully automated pipeline for the 
detection and the characterisation of microlensing-like 
flux variation s. We work in the "pixel-lensing" regime 
(|Gouldl 11991 . in which one looks for flux variations 
whose sources are not resolved objects, so that one has to 
monitor flux variations of every element of the image, fur- 
ther characterized by the fact that the noise is dominated 
by the underlying background level (the varying M31 sur- 
face brightness). As for this specific analysis, our strat- 
egy st arts from that described in ICalchi Novati et al.l 
(2005) with a few changes introduced to take into ac- 
count the peculiarities of the present data set. 

During the analysis we have to face two main sources 
of contamination: "fake" signals, namely spurious varia- 
tions to be attributed to cosmic rays, defects in the CCD, 
saturated pixels and so on; background intrinsically vari- 
able objects, that can either mimic microlensing signals 
or, somewhat more dangerously, add non-Gaussian noise 
to the light curves. 

Our pipeline can be schematically divided into four 
steps. First: detection of the potentially interesting flux 
variations. Second: characterisation of the light curve 
shape. Third: probe against the contamination by spuri- 
ous detections. Fourth: probe against the contamination 
by variable signals. 

3.1. Bump detection 

As for the first step we clo s ely fo l low th e strategy out- 
lined in ICalchi Novati et all (|2003l . l2005h . To begin, we 
detect flux variations along light curves using the L esti- 
mator (we ask L > 50 to get rid of too small S/N varia- 
tions). Each given flux variation enhances a signal over a 
few pixels (a "cluster") around the central one. We make 
use of the "Q" estimator to characterize the significance 
of the selected flux variations 11 . We fix a lower threshold 
Qth = 50. At this stage, therefore, we have to shift from 
the light curve analysis (the estimation of L and Q), to 
an analysis based on the spatial information across the 
CCD in which we have to distinguish, separate and pick 
up the flux variations associated to each different cluster. 
This search is somewhat biased in favour of light curves 
showing a single variation (in particular, all short period 
variables are in principle excluded). This first step is 
carried out using the (more numerous and less noisy) R 
band data only. 

3.2. Light curve shape 

The aim of the second step is to single out the varia- 
tions whose shape is compatible with that of a Paczyhski 
light curve. To this end, we use a series of selection 

11 It results L > whenever there are at least three consecutive 
points at least 3cr above the background. The value of Q is given 
by the ratio of the y 2 of a fla t baseline fit over that of a Paczyhski 
fit llCalchi Novati et al.|[200l) . 
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criteria. As a starting point we perform a 7-parameter 
Paczyhski fit in the two bands simultaneously 12 . As an 
output we retain the following parameters: the baseline 
levels; the time of maximum magnification, to; the full- 
width- at-half-maximum duration, ti/ 2 (proportional to 
the Einstein time, £e, multiplied by a function of the 
impact parameter Uq); the flux deviation at maximum 
(with respect to the baseline) , which we convert to mag- 
nitude and denote AR, together with the color of the 
variation, R — I, (AR is proportional to the source flux, 
<p* , multiplied by a function of the impact parameter); 
and the reduced \ 2 of the fit. We make use of the full 
parametrisation of the Paczyhski fit, looking for the Ein- 
stein time, the magnification and the unlensed source flux 
values, even if the intrinsic parameter degeneracy (linked 
to the fact that the source is not resolved), in most cases, 
does not allow one to accurately estimate the single pa- 
rameters. As a selection criterion we ask x 2 /dof < 10. 
This rather high threshold is motivated by the necessity 
to handle light curves contaminated by low level noise of 
non-Gaussian nature that can be attributed in particular 
to nearby blended intrinsic variables. 

As a second test on the shape, we ask for the bump 
to be suitably sampled by the observed data points. We 
split our analysis on the basis of a more or a less demand- 
ing requirement, so that we are going to refer to set "A" 
and set "B" of candidates, respectively. The details are 
given in the Appendix [AJ 

As a final test on the shape, we look at the charac- 
teristics of the detected flux variations. The two rele- 
vant parameters are the event duration, ii/2, and the 
flux deviation at maximum, AR. In order to appropri- 
ately delimit these parameter spaces, we have to balance 
for the efficiency, the expected event characteristics and 
the risk of contamination of the background of variable 
stars. We introduce a cut to exclude too faint variations, 
too noisy and therefore difficult to distinguish against the 
variable contamination. As a selection criterion, we ask 
for AR < 21.5. As for the duration, we do not expect mi- 
crolensing events to last more than 10-20 days fSect. ETTj) . 
Besides, we expect long-duration variations to be heavily 
contaminated by intrinsic variable signals. However, we 
prefer not to introduce any cut for this parameter allow- 
ing for long duration candidates. As detailed below, all 
of these are rejected in the following steps of the analysis 
anyway. 

3.3. PSF shape: spurious detection 

Up to now the analysis is based on the light curve 
pixel-photometry only (besides the initial "cluster" anal- 
ysis), in particular we do not make use of the PSF of 
the flux variations we are looking for. In this respect 
our approach is completely different from that based on 
the difference image analysis. As it is, however, this ap- 
proach suffers from a high risk of contamination by spu- 
rious variations. To reject them in an automated way, 
as a third step we perform an extremely rough difference 
image analysis. The underlying rationale is that, when- 
ever we detect a variation on a light curve, in order to 
retain it we want to "see" a corresponding well shaped 
PSF when we look at the image difference of the max- 

12 We use the CE RNLIB-MINUIT libraries, 
|bttp: / / cernlib . web . cern. ch/cernlib/ 1 



imum magnification minus the baseline. Futher details 
are given in the Appendix [B] 

3.4. Variable signals 

As a fourth and final step we probe the surviving vari- 
ations against the background of variable contamination. 
Our limited baseline, 50 days, does not allow us by it- 
self to carry out this programme. For this reason we 
make use of the 3-years baseline of the POINT- AGAPE 
data set 13 . The rationale is that we want to reject flux 
variations that show variability along the much longer 
INT baseline with a comparable flux deviation to that 
detected on our OAB data. As a first step, given the 
OAB detected variation, we look for the corresponding 
pixel within the POINT- AGAPE data set 14 . To probe 
variability along the INT light curve we use a Lomb pe- 
riodogram analysis. As an estimator, we use the powe r 
"Pr" as defined in Numerical Recipes (jPress et al.f l992). 
Second, we test whether the INT and the OAB flux devi- 
ations are compatible. To this end, given a relative flux 
calibration between the two data sets, first we rescale 
the OAB flux deviation; then we evaluate the difference 
of the observed flux deviations (the rescaled OAB minus 
the INT one), normalized by the INT error. We take this 
quantity, which we define 5(Af), as our estimator. As a 
selection criterion, we ask for Pr < 20 or S(Af) > 5. 

4. THE EXPECTED MICROLENSING SIGNAL 

In order to gain insight into the observed signal, first of 
all we have to evaluate the expected signal for our exper- 
iment al set up. We briefly out line our approach as fol- 
lows ([Calchi Novati et al.ll2005l) . First, we run a Monte 
Carlo simulation whose purpose is to give us the charac- 
teristics, and in particular the number, of the expected 
microlensing events. To this end, we have to specificy 
an astrophysical model and a model for the microlensing 
magnification, besides reproducing as closely as possible 
the actual experimental conditions. As for this last point, 
to account for those aspects of our pipeline that can not 
be accurately reproduced within the Monte Carlo, we 
carry out a simulation on the real data set of the events 
selected within the Monte Carlo. 

4.1. The model 

We consider the bulge and the disk populations of M31 
as sources, both M31 bulge and disk stars as lenses (we 
refer to these events as "M31 self lensing"), and MA- 
CHOs in both the M31 and the Milky Way dark matter 
halos. We assume an M31 distance of 770 kpc. 

As a model fo r the M31 luminous components we take 
the iKentl (|1989l ) bulge-disk decomposition, including the 
bulge ellipticity, and with the missing information of the 
vertical distribution of the disk modeled with a sech 2 law 
and scale height of 300 pc. For both halos we assume a 
spherical isothermal distribution with a core radius a = 
5 kpc. 

The peculiarity of M31 microlensing is that we look 
for flux variations of unresolved sources. The issue of 

13 Data collected at the 2.5m INT telescope du ring 1999-2001 
l|Auriere et al.|[200ll: IPaulin-Henriksson et lH]f2003T ) • 

14 The accuracy of the astrometric trasformation is below 1 pixel, 
but we must accept the limit given by the larger size of the OAB 
pixels, 0.58", with respect to the POINT- AGAPE ones, 0.33". 
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estimating the number of available sources is therefore 
particularly delicate. To this end, we consid er the value 
of the M31 surface brightness (as given by EES 1 1989) 
and the underl ying luminosity fu nction. As has already 
been remarked (|Ansari et a,l.ll 19971 ) . it turns out that only 
luminous sources (about Mj < 2) are expected to give 
rise to detectable events. In the most crowded region, 
this may sum up to hundreds of available source stars 
per pixel (this can be considered as a completely blended 
situation with respect to, for instance, studies towards 
the LMC). As for the luminosity function that we use 
to characterize the sources, it is worth stressing that we 
are most interested in its bright end (even if we need the 
information over the complete magnitude range for the 
normalisation), whereas, for the mass function that we 
need for the luminous lenses, we are rather interested 
into the opposite tail. 

We have made use of the IAC-star software 
(jAparicio fc Ga llart 2004) to build a synthetic luminos- 
ity functio n for the M31 bulge , following the procedure 
outlined in iBozza et alj (120081) . in particular as for th e 
metallicity distribution (|Saraiedini fc Jablonkal 120051 ). 
with the difference that we have now used, as a mass 
function, a power law £(m) oc mT a with index a = 1.33 
up to 1 M© and a — 2 above. For bulge lenses we have 
for consistency used the same mass function with upper 
bound fixed at on e solar mass. For the dis k luminos- 
ity function, as in iCalchi Novati et alJ ((2005), we make 
use of the local neighborhood dat a obtained by Hippar- 
cos corrected at the br ight end (|Perryman et alJ 119971 : 
Uahreifi fc Wielenlll997l ). For the disk lens mass function 
we follow as well the local determination (|Kroupall2007h 
with upper bound fixed at 10 M©. For MACHO masses 
we try a set of single values ranging from 10~ 3 to 1 M©. 

We fix th e total mass of the bulge to 4 x 10 10 M© 
(jKentl fl989| ) (this can be considered a "safe" value for 
microlensing analyses, for the purpose of evaluating the 
expec ted self-lensing sign al, as it is likely to be an upper 
limit ( Riff eser et al. 2006) f or this quantity), and that of 
the d isk to 3 x 10 10 M© (jKerins et al.ll200ttlRiffeser et alJ 
2006). Looking for microlensing effects, the value we are 
actually interested in is the stellar mass, and indeed our 
over all value for this qua ntity agrees well with the analy- 
sis of lTamm et al.l (|2007l ). We consider a uniform extinc- 
tion across the field, both foreground, E(B — V) = 0.062 
(|Schlegel et al.lll998D , and intrinsic (for this second term, 
this hypothesis should of course be taken only as a first 
order approximation), ext R = 0.19 for the bulge (jHanl 
fl99l iRiffeser et all 120061) and E(B - V) = 0.2 2 for 
the disk ((Stephens et al [ 120031 [Riffeser et all 12006). To- 
gether with t he M31 col or (not corrected for extinction) 
B-r = 1.3 (lKentlll989l) B - V = 1.0 and V - R = 0.8 
(|Walterbos fc Kennicuttl 119871 ). this translates into the 
values (corrected for extinction) M/Lr = 3.1, 1.1, for 
the bulge and disk respectively. 

The bulge velocity distribution is dominated by its dis- 
persive component, with line of sight velocity dispersion 
of d = 120 km s . For the disk we take into account 
both the dispersive motion, a = 60 km s^ 1 (a value 
that can be taken as un upper limit), and a circular 
bulk motion, with dis k circular velocity v = 250 km s~ 
(ICarignan et al.ll2006f). M31 proper mo tion is set accord- 
ing to Ivan der Marel fc Guhathakurtal (|2008h . 
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Fig. 1. — Characteristics of the expected events from the Monte 
Carlo analysis. From top to bottom: histograms for the distance 
from the M31 center, duration (full-width-at-half-maximum) and 
the flux deviation at maximum magnification. Histograms for event 
duration and flux deviation are given for self lensing only. For all 
of the histograms we show all of the events selected, namely for the 
full range of the parameters. In the bottom panel, the solid line 
indicates the threshold value we use in the pipeline. 

Given the values of the circular velocity, 220 km s _1 
and 250 km s _1 for the Milky Way and M31 respectively, 
we fix accordingly the central densities and therefore the 
overall masses of the halos and the values of the one 
dimensional velocity dispersions. As for the total dark 
matter halo mass value, within a truncation radius of 
100 kpc and 130 kpc for the Milky Way and M31 (we 
estimate the ratio of these values from those of the cir- 
cular velocities), we have respectively 1.0 x 10 12 M© 
(in good agreement w ith previous determination, e. g. 
IVallenari et al.l (|2006f ). but somewh at in excess wit h re- 
spect to the recent determination bv lXue et al.ll2008T ) and 
1.8 x 10 12 M©. 

4.2. The simulation 

Within the Monte Carlo simulation, given the astro- 
physical model, we generate microlensing signals with 
Paczyhski magnification corrected for finite size effects 
of the sources 15 , build the correspoding light curves and 
carry out a first rough selection, paying attention not to 
reject any light curve that could be detected through our 
pipeline. In particular, for a variation to be selected, as 
a unique criterion we ask for at least 3 consecutive points 
(with one point per night and reproducing the observed 
sampling) 3 sigma above the backrgound level (L > 
according to the estimator introduced in Sect. [3]). On 
the other hand, we are aware that we can not reproduce 
within the Monte Carlo, where we only deal with light 
curves, all of the the actual conditions of the pipeline 
we carry out in the real data set (where the analysis 
starts from the images). Amongst other effects, we most 
prominently can not reproduce crowding effects, the un- 
derlying variable signals and the sources of non-Gaussian 

15 For bulge sources we use the radius values from the IAC 
database, for disk sources we use a color temp erature relation eval- 
uated from the model of [Robin et al. (2003|) and we evaluate the 
radii from Stefan's law using a table of bolometric corrections from 
IMurdinl 120011) . For the micr olensing ampl i ficatio n we use the an- 
alytical expression derived in Witt & Mao (1994). 
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noise. Furthermore, we can not run the first essential 
step of bump "cluster" detection. To account for these 
effects we simulate those light curves that are selected 
within the Monte Carlo on the images (before the geo- 
metrical and all photometric corrections, namely, on the 
astronomical images just after the basic bias-flat fielding 
reductions) and then we run from scratch our pipeline. 
Conceptually, this is just a last step in the Monte Carlo 
that allows us to accurately evaluate the efficiency of our 
pipeline. The "efficiency", hereafter, should therefore 
be intended as that relative to the light curves selected 
within the Monte Carlo. 

Within the Monte Carlo each simulated event carries 
a (different) "weight" , Wi (where i is the index of the 
simulated events) that is linked in part to the drawing 
process and in part to the quantity we are evaluating (the 
microlensing rate), and is completely independent of the 
selection process . The expected number of events is 
therefore the sum of the weights, n oxp = Yli w ii where 
the sum runs over the events selected within the Monte 
Carlo (correspondingly we can estimate the associated 
statistical error based on Poisson statistics). Accord- 
ingly, by "efficiency" , e w , we mean the ratio of selected 
over simulated events, where the number we refer to is 
always given by the sum of the weights. This is usu- 
ally different (both numerically and from a logical point 
of view) from the actual ratio of selected over simulated 
events, a quantity that is not used even if in some cases 
it may be useful to be looked at. In particular, we do not 
expect, and in fact we do not find, these two values to 
be too different. Indeed, such a result should be taken as 
a hint of the presence of some bias in the way the event 
weights are distributed with respect to the selection pro- 
cess within the simulation. 

For each lens population, we simulate up to a few thou- 
sands events per field, with 500 events per field per sim- 
ulation in order to avoid overlap problems. Indeed, in 
particular in the inner M31 region, where we expect most 
of the events, simulated events may overlap (we draw at 
random from the Monte Carlo the events we simulate, 
with all their characteristics, including the line of sight) 
and thus lead us to bias the estimate of the efficiency 
of our pipeline. To test for this effect, for a fixed set 
of 500 events per field for which we had already evalu- 
ated the efficiency, we have carried out as many different 
simulations as needed, taking care to leave a minimum 
distance of at least 20 pixels among any couple of gener- 
ated events, so to exclude overlap problems. As a result, 
we have found no significant differences in the two anal- 
yses. Overall, we have simulated 12000 light curves to 
evaluate self-lensing efficiency, and 8000 for each value of 
the mass for MACHO lensing. 

According to the selection pipeline, in which we require 
for the flux variations to be large enough, AR < 21.5, out 
of the Monte Carlo we extract, and then simulate on the 
images, selected events with flux deviation at maximum 

16 As for this technical aspect, our analysis t herefo re differs, for 
instance, from that discussed in Kcrins ct al. (2001), as we draw 
all the values of the random variable according to their actual dis- 
tributions rather than according to the micro lensing rate ( f urther 
details are given in Appendix [C) . In addition, Kcrins ct alj (|2001h 
propose to perform the Monte Carlo simulation only to evaluate 
the pipeline efficiency whereas we make use of this tool also to 
evaluate the number of expected events. 
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Fig. 2. — Selection pipeline results: the sample shown here con- 
sists of the light curves selected after the PSF analysis. Scatter 
plot of the flux deviations versus the event duration (top panel); 
scatter plot of the Lomb "power" Pr versus the relative flux differ- 
ence (S(A/) (filled and empty symbols for set A and B variations, 
respectively). The squares indicate the two candidate microlensing 
events of our selection (according to the terminology introduced in 
Sect. 15.21 1. The solid lines delimit the excluded region (bottom 
panel, upper left corner) in the last step of the analsyis. The star 
symbol indicates the second event, beside Nl, belonging to both 
set A and set B. Bottom panel : for visualisation purposes we do 
not show a few rejected events with S(Af) < and Pr > 20. 

down to 17 AR — 21.8. This limit is used accordingly 
when we evaluate the number of expected events. This 
way we allow for the observed rms of the evaluated flux 
deviations versus the input values. The exact value of 
this threshold is not, however, essential, as long as we 
keep the selection criterion on the flux deviation fixed 
at AR < 21.5 coherently with the selection pipeline. A 
brighter threshold would translate into a larger value for 
the efficiency and, at the same time, a smaller number 
of expected events (not corrected for the efficiency), and 
vice versa. These two effects balance when we evalu- 
ate the number of expected events corrected for the effi- 
ciency. 

As for the expected characteristics of the observed 
events, in Fig. [T] we show the resulting distributions for 
the distance from the M31 center (both self lensing and 
MACHO lensing), and, for self lensing, the distribution 
for the durations (the full width at half maximum ti/2) 
and the flux deviations at maximum (AR). 

Finally, it is worth stressing that we simulate mi- 
crolensing events only. Therefore the simulation is re- 
stricted to saying whether and how our pipeline is go- 
ing to select microlensing signals but it can say nothing 
about whether it might select as a microlensing a variable 
signal of different origin. 

5. RESULTS 

In this Section we present and discuss the results of 
our analyses: the selection pipeline for microlensing flux 
variations and the evaluation of the expected microlens- 
ing signal. 

17 As a matter of fact, within the Monte Carlo, where we have 
only a statistical error, we find a much fainter "theoretical" lower 
bound for the flux deviation at maximum (Fig. [TJ. The efficiency 
analysis, on the other hand, showed us that the choosen treshold 
value for AR is appropriate, because the efficiency dramatically 
decreases when we consider too faint flux variations. 
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Fig. 3. — INT and OAB i?-band light curves extension of 4 
out of the 14 set A flux-variations OAB candidate microlensing 
events selected after the shape and the PSF analyses (left and right 
panels, respectively). The solid lines mark the evaluated baseline 
level, while the dotted lines show the maximum flux deviation value 
along the detected OAB flux variation (right panels) and this same 
value once rescaled on the INT data set (left panels) . For each INT 
data light curve we indicate the result of the Lomb periodogram 
analysis, the power Pjj, and the difference of the flux deviation 
along the OAB and the INT light curves, 6(Af). The units on the 
x axes are time in days (JD-2450000.0). The ordinate axes units 
are flux in ADU s~ x per superpixel. 



TABLE 1 
Microlensing selection pipeline 



THE RESULTS 





selection 


simulation 


criterion 


# events 


efficiency 


e w (%) 




set A set B 


set A 


set B 


bump detection 


4200 


41.2 ±3.5 


L x > 50 


3033 


35.5 ±3.2 


shape analysis 








X 2 /doi < 10 


2901 


32.7 ±3.0 


sampling 


174 241 


18.3 ±2.0 


21.9 ±2.4 


AR < 21.5 


75 108 


13.0 ± 1.5 


20.4 ±2.3 


PSF 


14 23 


13.0 ± 1.5 


15.7 ± 1.9 


variable 


1 2 


11.9 ± 1.4 


14.6 ± 1.8 



Note. — The results of the selection pipeline for mi- 
crolensing light curves: analysis and simulation. For each 
step we report the number of selected light curves and the 
efficiency of the pipeline (%) for the expected self-lensing 
signal. According to the choice of the sampling criterion, 
we have split our selection results in set A and set B (left 
and right column, respectively). 



5.1. The selection pipeline 

In TableQ]we report the results of the selection pipeline 
analysis together with the results for the efficiency of the 
corresponding analysis carried out on self-lensing simu- 
lated events. 

We start the analysis working over the complete set of 
pixels, namely 2 x (1340 x 1300) light curves. The ini- 
tial sample of selected flux variations consists of ~ 4000 
light curves. Within the shape analysis, the sampling cut 
severely reduces this initial set, and then the flux devia- 
tion cut leaves us with ~ 200 light curves, most of which, 
according to the PSF analysis, are to be attributed to 
spurious variations. Finally we are left with ~ 40 flux 
variations, divided into sets A and B (according to the 
sampling criterion), most of which we expect to be in- 



trinsic variables whose single-bump appearance is to be 
attributed to our short (50 days) baseline. (As outlined 
in Appendix set A flux variations are not a subsam- 
ple of set B: among those surviving the PSF cut, 14 and 
23 respectively, only two flux variations are in common 
between the two data sets, out of which one also sur- 
vives the last cut.) The results of the last cut analysis 
are shown in the bottom panel of Fig. [2] in the parameter 
space Pr-(5(A/). We find, as expected, the flux deviation 
of most OAB variations to be compatible with the cor- 
responding INT variations (small values of 5(Af)), for 
which, at the same time, we find a clear sign of variabil- 
ity (large value of Pr). As an example, in Fig.[3j we show 
the INT extension of 4 OAB selected light curves. Only 
for two selected OAB flux variations, instead, we find 
the corresponding INT extension to be flat. Therefore 
the selection pipeline finally leave us with two microlens- 
ing candidates (one belonging to set B only). The same 
sample of light curves is represented in the Ai? max — ti/% 
parameter space (top panel of Fig. [2|) . We find most of 
the variations located in the upper part (corresponding 
to small flux deviations), with set B variations (empty 
symbols) biased towards short durations. In particular 
we find the set A microlensing candidate (filled square 
symbol) located in a parameter space region where the 
contamination by the intrinsic variable signals is large. 
The only clear outlier is the bright and short set B mi- 
crolensing candidate. We discuss the selected candidates 
in detail in Sect. 15.21 

Comparing with the efficiency simulation analysis, a 
few points are worth being mentioned. First, it may 
look as if the "bump detection" step alone severely re- 
duces the overall efficiency. However, in fact, only a very 
small fraction of light curves that arc not selected at this 
point would pass all the other criteria. According to the 
same principle, the single cut that excludes most of the 
simulated light curves is the sampling criterion 18 . This 
is also the reason why we have split our selection pipeline 
at this level. As for the simulation, we stress that this 
simply reflects the choice we have made in the Monte 
Carlo to select light curves on the basis of the L > cri- 
terion only. Next, the PSF analysis proves to be a rather 
efficient criterion. Indeed it results that almost 50% of 
the simulated light curves fulfill this criterion and that 
this fraction rises to about 80% if we consider the sub- 
set of light curves that have already passed the bump 
detection and the shape analysis cuts. A usual reason 
that may cause the PSF Gaussian fit to fail is the pres- 
ence, near the simulated event, of some other resolved 
object. Often enough, however, in these cases also the 
pixel photometry we use may have problems. (For both 
of these related aspects, a proper difference image anal- 
ysis approach would be of course expected to give better 
results). Poorly sampled light curves, on the other hand, 
simply do not have enough points near maximum mag- 
nification, so that it turns out to be usually not possible 
to carry out a good enough PSF fit. Finally, as for the 
analysis on the POINT-AGAPE extension to check for 
variable signals, we have seen this cut to be essential to 
get rid of otherwise dangerous contaminating flux varia- 
tions. At the same time, this shows to be an extremely 

18 The overall efficiency would jump to ~ 20% taking into ac- 
count all the criteria except the sampling criterion. 
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TABLE 2 
Efficiency analysis: the results 
for macho lenses. 



mass (Mq) 


efficiency 


e w (%) 




set A 


set B 


1 


16.5 ± 1.2 


18.7 ±0.9 


0.5 


18.5 ± 1.1 


21.3 ± 1.2 


0.1 


16.5 ± 1.2 


19.9 ± 1.4 


10- 2 


13.8 ± 1.2 


16.5 ± 1.3 


10~ 3 


11.1 ± 1.3 


13.2 ± 1.4 



efficient criterion. Out of the complete set of simulated 
light curves, only about 13% of the INT light curves are 
clearly variable (Pr > 20) and this fraction falls to 7% 
when we add the demand for the INT flux deviation to 
be compatible with the OAB simulated one. 

In Table [2 we report the results for the efficiency of 
the simulations for MACHO lenses. With respect to self- 
lensing events there is a (rather small) effect linked to the 
different spatial distribution. The main effect is, how- 
ever, due to the value of the mass, which is linked to 
the event duration, with smaller values of the efficiency 
corresponding to decreasing values of the MACHO lens 
mass. 

In Tables Q] and [2] we have given the results for the 
efficiency as a single value for the overall set of simu- 
lated events. On the other hand, the efficiency does vary 
quite significantly for data binned, for instance, in the 
distance from the M31 center and/or in the flux devia- 
tion at maximum. We find the larger values for the effi- 
ciency, up to 30% or more, for bright events in the outer 
regions of M31. On the other hand, the expected number 
of events, not corrected for the efficiency, is larger near 
the M31 center for faint events (Fig. []}, namely, right 
where the efficiency is smaller (down to below 5%, de- 
pending on the choice of the binning). Overall, however, 
we find the expected number of events corrected for the 
efficiency to be rather insensitive, within the statistical 
error of the simulation, to any binning scheme, and this 
motivates our choice for the way to present our results. 

5.2. The microlensing candidate events 

The selection pipeline described in the previous section 
leaves us with two candidate microlensing events, which 
we name OAB-N1 and OAB-N2 ("N" indicates that they 
are both located in our "North" field) . Their characteris- 
tics are summarised in Tabled and their position within 
our field of view is shown in Fig. [If 9 . As for the sampling 
criterion, OAB-N1 fulfills both sets A and B demands, 
while OAB-N2 only the set B one. 

OAB-N1 : This is a relatively large flux variation, 
with significance bump estimators equal to L = 183 and 
Q = 108; quite short, ty% — 7.1 days, and not too bright, 
AR = 21.1. The OAB-N1 light curve is shown in Fig. OH 
together with its INT extension and the image difference 
around the candidate position upon which is based our 
PSF analysis. The INT data extension of the OAB-N1 
light curve appears to be flat (small Lomb periodogram 
power, Pr = 6.3, and significantly large value for the dif- 

19 The original M31 image has been taken from the CDS data 
base, http: / /cdsweb. u-strasbg.fr/ 




Fig. 4. — Projected on M31, we display the boundaries of the 
two 13'xl2.6' monitored fields. The filled circle mark the position 
of the selected microlensing candidate events. 



TABLE 3 

Characteristics of the microlensing candidates OAB-N1 
and OAB-N2. 



OAB-N1 OAB-N2 

OAB data source flux fixed 



a (J2000) 


Oh 42m 57s 




Oh 42m 50s 


8 (J2000) 


41°22'50" 




41° 


18'40" 


<^M3i (arcmin) 


7.1 






2.8 


t (JD-2450000.0) 


4433.8 ±0.2 


4467.3. 


-1.3 


4466.6 ±0.3 


h/2 (days) 


7 1+ 14 


2.6j 


-1.0 
-0.S 


2 56+ - 12 


AR 


21.1 ±0.2 


X9.1+0.S 


19.5 ±0.2 


R-I 


1.0 ±0.2 


1.1 ± 


0.1 


1.2 


X 2 /dof 


3.9 




1.4 


1.1 


t E (days) 




3.9 H 


-2.7 


4.181^1 


uo 




0.23-o.i2 




4>* R (ADU s" 1 ) 




15 


-12 


7.6 


<f>* T (ADU s" 1 ) 




28 


-22 


15.1 


R* 




20.2 H 


-1.4 


20.9 


I* 




19. 1 H 


-1.7 


19.7 



Note. — For OAB-N2, we show the results for the fit performed 
both with our data set alone (left column) and fixing t he source flux 
from a possible identification with a source in the IMassev et al.l 
(2006) catalogue. The error on the standard magnitude values and 
colors includes an extra ±0.1 mag term from the calibration equa- 
tion. 

ference of flux deviation, 6(Af) — 10.0). However, the 
quality of the Paczyhski fit is not good. A few points 
before the bump deviate significantly from the expected 
shape, and this is reflected in the rather poor value of x 2 - 
This might be attributed to underlying nearby variables, 
but we can not exclude it to be a sign of an intrinsic 
non-microlensing nature. Indeed, we have shown that 
OAB-N1 is located in a part of the t 1 / 2 — AR parame- 
ter space where the background of variable stars is large 
(top panel, Fig. A possible contamination, compati- 
ble with its "one bump" nature, besides a possible very 
long period variable, might come from some kind of erup- 
tive variable (even if Novae can be excluded because the 
flux variation is far too small). On the other hand, the 
descent is fairly well sampled and matches nicely enough 
the fitted Paczyriski shape. In the top right part of Fig. [5] 
we show the "color" light curve, namely the ratio R over 
/ band of the difference of the light curve flux, along the 
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Fig. 5. — The microlensing candidate event OAB-N1. Top 
left panel: R and /-band light curves (filled and empty symbols, 
respectively). The i?-band data are normalized with respect to 
/-band data using the color of the flux variation. The solid and 
the dashed lines are the best-fitting Paczyhski curves with to = 
4433.8 and t = 4430.6 (JD- 2450000.0), respectively (Sect. [531 for 
details). Top right panel: color light curve. Bottom panel (left): 
the extension along the INT data (the solid and dotted lines are as 
in Fig. [3J . The units on the x axes are time in days (JD- 2450000.0). 
The ordinate axes units are flux in ADU s _1 per superpixel (left 
panels) and magnitude (top right panel). Bottom panel right: OAB 
ij-band data surface plot image difference between the image at 
maximum magnification and the baseline level. The units on the 
space x — y axes is pixels, the surface plot values are in ADU per 
pixel. 

bump, and the background level, that we expect to be 
constant for microlensing. 

Through the pipeline, and in particular for OAB-N1, 
as an initial condition for the time of maximum magnifi- 
cation in the Paczyhski fit we choose the value of the time 
corresponding to the data point with the maximum flux 
value, in this case t = 4435.4 (JD-2450000.0). During the 
fit procedure it is difficult for this parameter to exit from 
the x 2 minimum well around this value (whose bounds 
are usually set by the sampling) and indeed as a result 
we find t = 4433.8 ± 0.2 (JD-2450000.0). Motivated by 
the OAB-N1 light curve appearance, irregular sampling 
and noisy aspect, we have therefore carried out a search 
for other x 2 minima in different regions of the to space. 
As a result we have found a new, well isolated, minimum 
with a lower value of the reduced x 2 (x 2 /dof = 3.5 to be 
compared with x 2 /dof = 3.9 found in the previous case) 
for t = 4430.6 ±0.4 (JD-2450000.0). Correspondingly, 
wc find new values for the duration and the flux deviation 
at maximum, tjy 2 — 18 days and AR = 21.5 (this last 
value being at the limit of our threshold cut), namely, 
a rather different result from the previous one. Such a 
light curve would have still been selected as a microlens- 
ing candidate within our pipeline. However, the longer 
duration would have futhcr weakened its microlensing 
interpretation (Sect. [5~4|) . The co-existence of these two 
minima might be suggestive, accepting the microlensing 
hypothesis, of a binary lens or binary source solution. 
The available data, however, do not allow us to robustly 
test this hypothesis 20 . 

20 Through the selection pipeline we have only considered (sin- 
gle bump) Paczynski like microlensing variations. The analysis of 
the possible binary lens solutions for OAB-N1, together with a sys- 
tematic search for binary-like flux variations, will be presented in 



4462 4464 
residuals 



■in 

4420 



4430 4440 4450 44 

OAB R&I light curve 



4460 4462 4464 4466 

OAB R-I light curve 



50 , I 

m y. *n 

2 5 L , \ , , 

1400 1600 1800 2000 2200 
INT R light curve 



Fig. 6. — The microlensing candidate event OAB-N2. Top panel: 
normalized R and /-band light curves (filled and empty symbols, 
respectively). The ft-band data are normalized with respect to I- 
band data using the color of the flux variation. The solid line is 
the best-fitting Paczynski curve. In the insert we show the residual 
with respect to the Paczynski fit along the bump. Bottom panels, 
left: color light curve; right: extension along the INT data (the 
solid and dotted lines are as in Fig. EJl. The units on the x axes are 
time in days (JD-2450000.0). The ordinate axes units are magni- 
tude for the bottom left panel, and flux in ADU s _1 per superpixel 
for the remaining panels. 

OAB-N2 : This looks (Fig. EJJ) like an extremely bright 
and short flux variation (ijya = 2-6 days, AR — 19.1), lo- 
cated near the M31 center (at a distance d = 2.8'), with 
a completely flat INT extension (P R = 8, 5{Af) = 25). 
The sampling along the bump is, however, extremely 
poor (indeed, both the time of maximum amplification 
and the flux deviation at maximum are not strongly con- 
strained, Table [3]). In particular, the lack of data points 
in the descent prevents us from probing the expected mi- 
crolensing symmetric shape, so that it is difficult to draw 
definitive conclusions about its nature. On the other 
hand, the rising part of OAB-N2 looks achromatic, and 
is rather well constrained, and this strengthens the mi- 
crolensing hypothesis. Indeed, both the shape and the 
achromaticity are hardly compatible with the only pos- 
sible kind of contamination for such a kind of flux vari- 
ation, if not microlensing, namely some sort of eruptive- 
like object. Furthermore, though few, the points along 
the bump allow us to get to a reasonable fit for all the 
parameters of the microlensing event. At the 1 a level, 
the x 2 analysis gives us best values and lower bounds for 
the impact parameter and the source flux and best value 
and upper bound for the Einstein time (with, however, 
a rather large relative 1 a error, even exceeding 50%, 
Table [3]). Finally, the guess for the flux source value 
allows us to carry out a more constrained color analysis 
(as for OAB-N1, but subtracting the source flux from the 
baseline, bottom left panel of Fig. [6J. As for the value 
of the source flux, the results of the Paczyhski fit have 
been confirmed by the cross-ide ntification of t he po ssible 
source star in the catalogue of iMassev et al.l (|2006f ). In- 
deed, within 1 pixel of our evaluated position, we find a 
typical red giant with Rc — 20.9 and R — I = 1.2 (val- 
ues fully compatible with those evaluated from our data 
set alone). The knowledge of the source flux allows us 

IBozza et all pOOlf) . 
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TABLE 4 

The expected number of microlensing 

EVENTS. 





n 


cxp 




set A 


bUT, id 


bulge-bulge 


0.35 ±0.07 


0.44 ± 0.08 


bulge-disk 


0.18 ±0.02 


0.23 ± 0.03 


disk-bulge 


0.06 ±0.01 


0.07 ±0.01 


disk-disk 


0.015 ±0.002 


0.019 ± 0.002 


SL 


0.6 ±0.1 


0.8 ±0.1 


mass (Mq) 






1 


0.7±0.1 


0.9 ±0.2 


0.5 


1.1 ±0.2 


1.3 ±0.2 


0.1 


1.6 ±0.1 


1.9 ±0.1 


lO" 2 


2.0 ±0.2 


2.3 ±0.2 


10~ 3 


1.3 ±0.2 


1.5 ±0.2 



Note. — The expected number of mi- 
crolensing events, for M31 self lensing (lenses 
belonging to either the M31 bulge or disk) 
and MACHO lensing (lenses belonging to ei- 
ther the M31 or the Milky Way halo), for 
a full halo with different values of the MA- 
CHO mass. For self-lensing events we report 
also the number for each lens-source popu- 
lation we consider. We report the results of 
the Monte Carlo simulation corrected for the 
efficiency evaluated according to both set A 
and B criteria. 

to better constrain the fit parameters, in particular the 
time of maximum amplification and the flux deviation at 
maximum (in fact t$ shifts back to the position of the 
last observed data-point so that, accordingly, AR gets 
fainter). Furthermore, we may now completely break 
the degeneracy among the amplification parameters. In- 
deed, we estimate t-g = 4.2+q; 4 days and u a = 0. 22^0^3 

(^4max = 4.5). 

We have also searched for X-ray counterparts in 
the XMM-Newton archive 21 and verified that neither 
of the microlensing candidates' coordinates correspond 
to any of the ident i fied sources in the EPIC images 
(IStriider et alj 120011 : iTurner et all 120011: IShirev et all 
1 2 lh . In the case of OAB-N2 there is an X-ray counter- 
part, but only within ~ 30", so that we can safely rule 
out the identification. 

5.3. The expected number of events 

The Monte Carlo simulation, completed by the effi- 
ciency analysis, allows us to estimate the expected num- 
ber of microlensing events for our experimental setup. In 
Table|4]we report the resulting values, already corrected 
for the efficiency. The results we obtain using the set B 
criteria is larger by up to about 30% with respect to set 
A. This is a combined effect of the larger value of the 
efficiency and of the effective baseline length increase. 

As for self-lensing events, it turns out that most of 
them, ~ 50%, are due to bulge-bulge events (lens- 
source, respectively), with the remaining distributed al- 
most equally between disk-bulge and bulge-disk events. 
The second configuration is enhanced in the South field 
because in this case we see the bulge in front of the disk. 
It is also worth noting that, according to our model, al- 
most half of the overall bulge mass is located within our 
cone of view, but only about 1 /7 of that of the M31 disk. 

21 http://xmm.esac.esa.int/xsa/. 



As for MACHO lensing, about 2/3 of the events are 
expected to belong to the M31 halo. Actually, the overall 
mass within our cone of view of the M31 halo is larger 
than the Galactic one by a factor of a few thousands, but 
Milky Way halo lensing is strongly enhanced by the much 
larger size of the Einstein radius. We expect the largest 
number of MACHO lensing events from 10 _1 — 10 -2 M Q 
objects, as our efficiency drammatically drops for smaller 
MACHO masses. 

Together with the expected event number we report 
the statistical error associated with the simulation on the 
basis of Poisson statistics. The uncertainties intrinsic in 
the model (whose detailed discussion is beyond the scope 
of the present paper) make, however, the corresponding 
systematic error in principle even larger. 

5.4. Discussion 

We have designed our (fully automated) pipeline so as 
to get rid of most forms of stellar variability, while pre- 
serving genuine microlensing. As a result, therefore, all 
survivors could in principle be variables but with vary- 
ing (hopefully large) degrees of confidence that they are 
not. Coming to the present data set, we have selected 
two flux variations compatible with a microlensing sig- 
nal. While discussing the possible physical meaning of 
this result, however, we have to keep in mind that, al- 
though with different degrees of confidence, we have no 
strong evidence in favour of the microlensing hypothesis 
for either candidate: OAB-N1 because of both the light 
curve appearance (reflected in the x 2 value) and its posi- 
tion in the duration/flux-deviation parameter space, and 
OAB-N2 because of its extremely poor sampling (that 
has not prevented us, however, from finding a convinc- 
ing microlensing Paczyhski fit further confirmed by the 
possible identification of the source) . In conclusion, look- 
ing at the candidate light curves, we might be tempted to 
classify OAB-N1 as a very "poor" candidate, and OAB- 
N2 as a "good" one, keeping in mind, however, that this 
might simply reflect the bias induced by the observa- 
tional sampling. 

With the care suggested by the above discussion, we 
come now to the comparison of the observed events with 
the expected signal. As for the event characteristics we 
may take as a reference the distributions shown in Fig. [TJ 
Given the caveat that these distributions can not be com- 
pared directly with the results of the analysis because of 
the correlation of the pipeline efficiency with the event 
characteristics, both candidates appear compatible with 
the expected signal with respect to the duration and the 
distance from the M31 center (being both shorter and 
nearer to the M31 center, OAB-N2 appears, also in this 
respect, a stronger candidate than OAB-N1). Coming 
to the number of events, we must compare our observed 
candidates with the results reported in Table [4j First, 
we are bound to consider the case that MACHO lensing 
does not contribute at all, so that we are left with the 
expected lensing signal from the luminous populations 
alone, with an expected number of events for self-lensing 
somewhat smaller than 1. Even allowing for both our 
candidates to be genuine microlensing, this is still, ac- 
cording to Poisson statistic, fully compatible with the 
observations. 

On the other hand, the expected number of MACHO 
lensing events is not large if compared to that of self 
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lensing. Even for a full MACHO halo we would expect 
about as many self-lensing events as MACHO lensing. 
It is clear, therefore, that on the basis of the number 
of events alone our statistics are still too small to draw 
conclusions on this contribution. 

It may be asked, then, how we might disentangle the 
two signals: self lensing from MACHO lensing. First, 
an improvement in the event statistics might help us to 
further exploit the expected differences in the spatial dis- 
tributions (top panels, Fig. [TJ). Second, a good enough 
sampling (such as was not the case for the present se- 
lection) might allow a sufficient characterisation of the 
selected flux variations. In turn, this could give possible 
hints on the nature of the lens, as was the case for the de- 
tailed analysis linked to the finite source siz e effect for the 
PA-S3 event carried out by iRiffeser et alj (|2008t ). Both 
these reasons motivate the need for further observations. 

6. CONCLUSIONS 

In this paper we have presented a thorough analysis 
of the data acquired during the second season of our 
pixel-lensing observational campaign carried out towards 
M31 with the 152 cm Cassini telescope located in Loiano. 
We have developed a fully automated pipeline analysis 
for the detection and the characterisation of Paczyiiski 
like microlensing flux variations. We have evaluated, 
through a Monte Carlo simulation, the expected signal, 
given an astrophysical model for M31 and the efficiency 
of our pipeline. As a result we have selected 1-2 candi- 
date microlensing events. This output turns out to be 
in fair agreement with our evaluation of the expected 
self-lensing signal from the luminous M31 components. 
As for the would-be MACHO lensing signal, the statis- 
tics of events are still too small to draw firm conclusions. 
Further observations might help us to disentangle be- 



tween the two populations. First, because of the larger 
statistics of events, second (provided a better sampling) 
because of a better characterisation of the light curves. 

Previous analyses, in particular those carried out us- 
ing data collected at the INT telescope by the POINT- 
AGAPE and the MEGA collaborations, have reached dif- 
ferent conclusions on the dark matter halo content in the 
form of MACHOs a l ong the line of sight towa rds M31 
(jCalchi Novati et all 120051 : Ide Jong et all l2006h . Their 
fundamental point of disagreement may be traced back 
to the evaluation of the expected self-lensing signal. In 
this perspective our campaign, as well as that carried 
out by the A NGSTROM collabo ration right towards the 
M31 center ()Kerins et alj [20061. may help in shedding 
more light onto this important issue. 
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APPENDIX 

THE MICROLENSING PIPELINE: THE SAMPLING CRITERION 

The rationale of the present cut is to test whether the detected flux variation is sufficiently sampled. As a starting 
point we use the results of the Paczyhski fit, and in particular the values of the time at maximum magnification, to, and 
of the bump duration, the full- width- at-half- maximum ti/2- Given these values we identify six intervals along the light 
curve, symmetric around to, [t<t -3 tx/2], [to - 3 t 1/2 < t < to -ii/2/2], [t -t 1/2 /2 < t < to], [to < t < to+h/2/2], 
[t +t 1/2 /2<t <t + 3t 1/2 ] and [t > t + 3 t 1/2 }. 

As a first set of criteria we ask for at least n m ; n data-points in at least 3 out of the 4 inner intervals, within to ±3 ti/2, 
and at least 1 data-point in both the joined outer intervals t < to — fi/2/2 and t > to + £1/2/2, namely, we allow the 
same data-point to be counted twice. The value of n m i n is fixed according to the duration, with n m i n = 1,2,3 for 
ti/2 < 5, < 15, > 15 days, respectively. 

As a second set of criteria we ask for at least 1 data-point in at least 2 out of the 4 inner intervals and at least 1 
data-point in at least 1 of the two tail intervals (t < t — 3 t 1 / 2 and t > to + 3 £1/2)- 

We refer to the two set of criteria, and to the corresponding selected flux variations, as set "A" and set "B", 
respectively. 

The set A criterion is more demanding as for the sampling along the inner part of the flux variation, whereas 
set B makes, in a way, a stronger demand on the coverage of the far tails, so that set A selected events are not a 
simple subsample of set B events. Furthermore, it results that the time of maximum magnification for set A events is 
constrained within the limits of the observational run whereas for set B it can also fall (slightly) outside. Set B flux 
variations therefore enjoy an effective longer baseline. 

The sampling analysis is carried out along i?-band data only. 

THE MICROLENSING PIPELINE: THE PSF CRITERION 

As outlined in Sect. [3) through this criterion we want to check whether the detected variation along the light curve 
can be attributed to a physically meaningful flux variation (a variable star or microlensing) or to some kind of spurious 
signal (cosmic rays, bad pixel, seeing effects and so on). To this purpose we have to turn to an analysis carried out 
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on the images. As a criterion we test whether the spatial form of the detected bump has a well enough shaped PSF. 
As we are not carrying out difference image photometry, we can safely stop at the first order approximation of a 
2-dimensional Gaussian, namely 



P(x, y) = exp ■ 



1 



2(1 -P 2 ) 



Vx\ , (V-Hy\ _ 2p ( x - fi x \ ( y - n v 



(Bl) 



We carry out this analysis on a difference image (a window of 19x19 pixels around the central pixel Xo,ya i n which 
the flux variation has been detected) built as follows. Given the time of maximum magnification of the bump variation, 
to, we look for the nearest night of observation, and we average the images collected during this night, taking note of 
the average seeing of the night. As a baseline image we take the average of a set of images chosen along the baseline 
with a similar seeing. 

As a criterion, we ask for the fit with respect to the 2-dimensional Gaussian to converge "properly" 22 . Besides the 
actual demand for the fit to converge we also ask the following : \xq — (i x \ < 1.5, \yo — H v \ < 1.5, 1 < o x < 3, 1 < a y < 3 
and \p\ < 0.5. This additional set of criteria is necessary as it is very easy, in case of spurious variations, for the fit to 
converge while it slips at the border of the physically meaningful bounds. 

It is worth noticing that a different parametrisation for the 2-d Gaussian profile with respect to that chosen in Eq. IB 1 1 
is sometime used, namely 



p l(x,y) = exp 



(x cos <fi + y sm 4>) 2 (x sin 4> — y cos < 



2*2 2a 2 



(B2) 



In fact the two formulations are equivalent, one can easily evaluate the parameters cr a , ah and as a function of a x , <r y 
and p, and vice versa (it must be noted, however, that the parameter p plays a double role, giving both an inclination 
and a distortion). Within the parametrisation of Eq. IB2I we find our cuts to roughly correspond to 1 < a a < 3, 
1 < <7& < 3, and <j> w arbitrary. Still, we prefer Eq. IBll as the fit procedure looks, in this case, more stable. 



MONTE CARLO SIMULATION: THE DRAWING PROCESS 

In this Appendix we discuss some further details regarding the Monte Carlo simulation described in Sect.Q] 
A microlensing event is fully specified by 10 parameters. The line of sight, specified by two angles, 9, p, out of which 
we determine the (angular) position within our fields as x = pcos(#),y = psm(d); the lens distance D\ and mass pn\ 
the source distance D s and flux <p s ; the lens-source relative velocity, specified by a modulus v r and a phase ip v ; finally 
the microlensing amplification parameters, the impact parameter uq and the time of maximum amplification to- Each 
of these parameters has its own physical distribution function from which we draw a value for each event realisation. 
Even if we simulate only R light curves, the available information within the luminosity functions used allows us to 
evaluate also the source radius, so to properly take into account the finite size source effect, as for the microlensing 
amplification. 

A key aspect of our Monte Carlo scheme is the "weight", Wi, that we associate to each event. This is linked in part 
to the drawing process and in part to the physical process we are studying. As for the first part, this follows from 
the fact that for a given distribution function we have two choices as for the drawing process. Either we may directly 
draw according to the distribution, and in this case the corresponding weight is Wi — 1, either we may draw with a 
uniform distribution and give the event a weight proportional to the distribution function. As for the second part, we 
take the weight to be proportional to the following 

Wi oc (2 R E Umax v r AT bs) x n s , (CI) 

where i?E is the Einstein radius, Umax the maximum value for the impact parameter, AXobs the duration of the 
observational campaign, n s the number surface density of available sources (fox a fixed line of sight). In fact, according 
to our simulation scheme, first we fix the lens line of sight and specify all of the events characteritics. Then , as the 
purpose of the Monte Carlo is to evaluate the number of expected microlensing events, we have to evaluate how many 
sources (given the M31 surface brightness and luminosity function) are available within a surface delimited by the 
event characteristics and the experimental conditions. Indeed, the term 2 i?E Umax v r AXobs m Eq. ICll is the surface 
area on the lens plane (for a fixed event configuration) where it is possible to find a suitable source. 

The last contribution to the weight comes from the need to evaluate the number of available lenses. For each lens 
population this is computed as the ratio of the total mass of the component within the observed field of view divided 
by the average lens mass. More specifically, the weight linked to the choice of the line of sight is normalized to the 
total mass. 

As for the drawing of the line of sight we have the following: for lenses in the Milky Way halo we assume the 
spatial distribution across the observed field of view to be uniform; for lenses in the M31 halo, because of the assumed 
spher i cal sy mmetry, the resulting distribution for the angl e 9 is uniform ; for bulge and disk lenses we make use of the 
iKentJ (|1989f ) bulge-disk decomposition (whenever using the [Hi] (fl989h profile we use a uniform distribution and then 
attribute to the event a "weight" as described above). 

Once the event configuration, and its weight, have been specified, we build the corresponding light curve and, 
given our observational set up, we evaluate whether the event is selected or not. This selection process is completely 

22 We make use of the HFITH function of the CERNLIB. 
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TABLE 5 

Monte Carlo parameters: distributions function and limits for each lens/source population. 





MW halo 


M31 halo 


M31 bulge 


M31 disk 


limit 


e 


uniform 


uniform 


K89 


K89+/«) 


field of view 


p 


uniform 


PHM3l(*\) 


K89 


K89+/(C) 


field of view 


A 




PHM3i(r) 


K89 


K89+/(C) 


model 




SM 




It-" 


It-" 


model 


D a 






K89 


K89+/(C) 


model 


4>s 






IAC 


Hipparcos 


Mj < 2 


V r 


Gaussian 


Gaussian 


Gaussian 


Gaussian 


0,oo 




uniform 


uniform 


uniform 


uniform 


0,2tt 


"0 


uniform 


uniform 


uniform 


uniform 


0,«MAx(= 2) 


to 


uniform 


uniform 


uniform 


uniform 


0, AT obs 



Note. — The parameters have been introduced in Appendix \C\ The details for the den sity function of the 
Milky Way and the M31 halos and for the source luminosity functions are given in Sect. 14.11 K89 stands for the 
Kent ( 1989) bulge-disk decomposition. For the disk, as only a two-dimensional distribution on the plane is given, we 
consider a vertical distribution /(C) = cos h(— (C/Cft)) 2 ; where is the vertical scale height. The index for the power 
law mass function are given in Sect. 14.11 In the last column we report the range of v ariability for each parameter: 
for 0,p "field of view" indicates the two observed fields, "Nord" and "Sud" (Sect. |2~TT |; for the distances we take the 
physica l bou ndaries of the component considered, in particular for the halos we use a truncation radius as specified 
in Sect. |4~T1 for the bulge ~ 7 kpc around the M31 distance while for the disk vertical component we draw directly 
from the given distribution (for a fixed line of sight this allows us to evaluate the physical distances); for bulge and 
disk lenses fi t is drawn with a lower limit of 0.08 Mq and upper limit of 1.0 Mq and 10.0 Mq, respectively. 

independent of the event weight. The selected events delimit the overall integration space, so that the expected number 
of events is given by the sum of the weights of these selected events only. (This same set of events is then taken, with 
all of their characteristics, as the input set for the simulation on the images described in Sect. 14.21 ) 
In Table [5] we report, for each parameter, the distribution function used and its range of variability. 
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